Restoration of kx factorization for low px hadron hadroproduction 
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We discuss the applicability of the kr factorization theorem to low-pr hadron production in 
hadron-hadron collision. It has been shown that the kr factorization for high-pr hadron hadropro- 
duction is broken by soft gluons in the Glauber region, which are exchanged among a transverse- 
momentum-dependent (TMD) parton density and other subprocesses of the collision. We explain 
that the contour of a soft loop momentum can be deformed away from the Glauber region at low 
Pt, so the above soft gluons are factorized by means of the standard eikonal approximation. The 
kr factorization is then restored in the sense that a TMD parton density maintains its universal- 
ity. Because the resultant soft factor is independent of hadron flavors, experimental constraints 
£SJ ■ on its behavior are possible. The kr factorization can not be restored for the transverse single- 

spin asymmetry in hadron-hadron collision at low Pt, since the leading-order soft divergence is not 
factorizable. 
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The kr factorization theorem has been widely applied to inclusive and exclusive processes in perturbative QCD 
PHI]. This theorem holds for simple processes, such as deeply inelastic scattering (DIS) and Drell-Yan production. 
Recently, it was found that the kx factorization breaks down for complicated processes like high-p^ hadron production 
in hadron-hadron collision 0-Q 

0) . Hi(pi) + H 2 ( P2 ) -> H 3 ( P3 ) + Hi{ Pi ) + X, (1) 

. 

where pi denotes the momentum of the hadron Hi. The kinematic region, in which the produced hadrons are almost 
t^j- , back-to-back azimuthally (relative to the collision axis) , was analyzed in a model field theory containing scalar particles 
J> 1 and gluons under an abelian gauge group. The kx factorization, if applicable, gives the leading-power differential 
cross section [7H9( 

in . 

E ^ Ei d 3 p 3 d 3 p 4 = ^ J d(J i+]^k+if l /ifj/2dz/kdi/i, (2) 

which involves the transverse- momentum-dependent (TMD) parton densities fi/H [loLITlj . the fragmentation functions 
dff n, and the parton- level differential cross section dcri+j^k+i ■ The sum over the flavors and the integral over momenta 
of the partons are implicit in the above expression. 

When factorizing the TMD parton density fin, infrared divergences from the soft gluons were identified, which are 
exchanged among fan and other subprocesses of the collision [l,|9j]. These divergences, violating the universality of 
fi/i, break the kx factorization for the hadron hadroproduction. The source of the factorization breakdown is briefly 
explained below. Consider the one- loop diagrams in Figs. QJa), [T{b), andQJc) with radiative gluons being emitted by 
the spectator in Hi and attaching to the active partons in H 2 and H 3 , respectively. These active parton lines 
can be eikonalized, if focusing on the collection of infrared divergences. The eikonal line from H3 corresponds to the 
Wilson line appearing in the operator definition of /j/i- The infrared divergences in Figs.QJa) andQJb), which carry 
the information of the TMD parton density fj/2, need to cancel in order to ensure the universality of fi/\. However, 
the eikonal propagators associated with the partons in Hi and H4 are summed into an imaginary piece 

= -27ri<5(Z+), (3) 



-l+ + ie l++ie 

where the first (second) term comes from Fig. [UJa) (HJb)), and I denotes the loop momentum. The 5-function leads 
to a residual soft divergence from the Glauber region with l + = [l2l |. 
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(a) (b) (c) 

FIG. 1: Some one-loop diagrams relevant to the factorization of /,yi in hadron hadroproduction. Hermitian conjugates of these 
graphs also contribute. Particles in the above diagrams do not carry colors in an abelian gauge group. In the unpolarized case, 
both the solid and dashed lines represent scalars. In the case of single-spin symmetry, the solid lines represent fermions and 
the dashed lines represent scalars. 



This one- loop residual soft divergence may not cause trouble, if the leading-order (LO) amplitude is real. The 

expansion of the differential cross section for the hadron hadroproduction up to next-to-leading order (NLO) gives 

\M\ 2 = \M (0) \ 2 + 2Re[M (0) M w *}. (4) 

According to Eq. ([3]), the residual soft divergence will be purely imaginary, if is real, so it does not survive in 

the second term RefA'l^.M^*]. At two loops, the residual soft divergence, arising from the real product S(l t)S(lt) , 
does exist in the kr factorization. Certainly, the product <5(Z^")<J(Z^") does not break the collinear factorization |KML6l |: 
it has been explicitly demonstrated Q that the two-loop residual soft divergence cancels, when the parton transverse 
momenta are integrated out. The cancellation occurs between the diagrams with the two soft gluons on the different 
sides of the final-state cut and the diagrams with the two soft gluons on the same side. That is, the universality of 
a parton distribution function defined for the collinear factorization is maintained in the hadron hadroproduction. It 
was also found that the Glauber effect exists in complicated dijet production in hadron- hadron collision [17] . but does 
not in simple Drell- Yan processes [l8[ . 

In the present work we shall have a closer look at the failure of the kr factorization for the hadron hadroproduction. 
If the residual soft divergences could be factorized from the collision, the universality of fan would be recovered at 
the price that Eq. ([2]) involves an additional soft factor. However, the soft gluons are characterized by momenta in 
the Glauber region as stated before, in which the spectator line of Hi can not be eikonalized 0, 0. Therefore, the 
Ward identity argument used in standard factorization proofs [l9|, [2(| does not apply. It will be shown that the kr 
factorization actually holds at low pr (e.g., few GeV at Tevatron), though it breaks down at high pr as claimed in 
Q . The argument is that the contour of a soft loop momentum can be deformed away from the Glauber region at low 
pr, so the eikonal approximation applies to the spectator line of H±. The soft gluons responsible for the factorization 
breakdown are then factorized into a soft factor, and the kr factorization is restored in the sense that fy\ maintains 
its universality. Since the soft factor is independent of hadron flavors, experimental constraints on its behavior from 
some processes are possible, based on which predictions for others can be made. 

At low pr, the dominant contribution to the hadron hadroproduction comes from the region with a small parton 
momentum fraction xi = kf Jp\ <C 1, in which kyr is not negligible in the hard scattering either (2lT |. Another process 
to which the kr factorization applies is the semi-inclusive DIS [22]. In this case the small x region is reached by lowering 
the square of the transverse momentum of the outgoing hadron with respect to the virtual photon direction. The 
residual soft divergence from Figs, [lja) and tUb) is collected by 

T W = 2Afl 4 /— 2(p+ -k+)(k 2 + k i )-(ki+k 3 )6(l+) 

L 5 J (27r) 4 (fc 1 -0 2 (pi-fci+0 2 ' 2 (-2fc 1 + fc 3 --|k 1T -k3T-lr| 2 ) ' [ ' 

where only the piece relevant to the soft factorization is shown explicitly. A and g denote the triple-scalar coupling 
and the gluon-scalar coupling, respectively [ill ]. The transverse loop momentum lr in the numerator, being smaller 

than other terms, has been dropped [2l|, [23j|. The Glauber region involved in Tjp is defined by = lr) ~ 

(0,A 2 /E, A) with the center-of-mass energy E and the small hadronic scale A ~ k\r- Because of I 2 = — l T , the 
diagrams with real gluon emissions were not included in Fig. [T] [9] . 
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In the Glauber configuration the spectator propagator l/(pi — fei + I) 2 can not be eikonalized into l/[2(pi — kx) ■ I], 
since the term from the transverse momentum kiT — 1t| 2 ~ A 2 , being of the same order as (p^ — k^)l~ ~ A 2 , should 
be kept. To make our statement explicit, we consider the two poles in the l~~ complex plane for Eq. ([5]), 

r- h - |!T-k 1T | 2 |iT-k 1T | 2 

' ~ 1 2fc+ ' " fcl + 2(p+-fc+) (6) 

with k± ~ A 2 /E being demanded by the on-shell condition (pi — fci) 2 = 0. For a large xi, i.e., for fc^ ~ , these 
two poles in different half planes are close to each other, such that the contour of l~ is pinched and must go through 
the region of l~ ~ A 2 /E as shown in Fig[2][a). Therefore, |kiT — 1t| 2 and (pf — ki)l~ are of the same order, and the 
cikonal approximation fails as stated above. Picking up the latter pole in Eq. ([B]), Eq. (J5J) leads to 

T (1) - -A 2 -\ 4 f d2lT 1 ~ Xl (fe 2 + fc 4 ) ■ (fci + fc 3 ) 

47T iAg y ^ |kiT _ ^ _ 2fc + fc _ _ _ k3T _ lr|2 • • • ■ l<) 

For a small x\ or fc^ <C pf , the first pole moves away form the origin and becomes located at Z - ~ A 2 /fc 1 f , while the 
second one remains of 0(A 2 /E). One can then deform the contour of Z~ in the complex plane, so that the region of 
l~ ~ A 2 /S is avoided as shown in Fig. [2jb). We then have the hierarchy (pf — k^)l~ ~ (p^/fc^A 2 » |kix — 1t| 2 ~ A 2 , 
and the eikonal approximation can be justified for the spectator line in H\. To verify our argument, we compare the 
result from Eq. (O and that from the simplified integral with the eikonal approximation, 2(p^ — k^ )/(pi — k\ + Z) 2 sa 
1/Z-, 

r (i)«fe = 2 , 4 [jft (fcg + fcj) ■ (fci + fc 3 )^+) 

L ~ ff 7 (2 7 r)4(fc 1 -Z) 2 Z-Z 2 (-2fc+fc 3 --|k 1T -k3T-lT| 2 ) ' [ ' 

The pole Z~ = — ie from the eikonal propagator, corresponding to the second pole in Eq. (|6]), gives 

T (i)eik = _. 2- A 4 f d2 h 1 (fc 2 + fc 4 ) • (fci + fc 3 ) (q) 

L 9 J (2tt)4 (|k 1T - 1 T | 2 - 2k+K)l 2 T -2k+k~ - |k 1T - k 3T - 1 T | 2 " ' ' { ' 

where 2&i fcf can be reexpressed as x\k\ T j(x\ — f) after employing (pi — fci) 2 = 0. It is obvious that Eqs. (O and 
([9]) are quite different in the large X\ ~ O(l) region, but identical, if the terms proportional to x\ are negligible. 

For real QCD processes, the spectator lines in Fig. [T] are replaced by gluon lines. There could be infinitely many 
spectator lines, ie., rung gluons, forming the so-called ladder diagrams. The diagrams with the Glauber gluons 
attaching to rung gluons generate the residual soft divergences. It has been known that the region with strong 
rapidity ordering for the ladder diagrams give s a dominant contribution, which has been summed into the Balitsky- 
Fadin-Kuraev-Lipatov evolution equation [24j j . The strong rapidity ordering corresponds to ki <C Pi in the model 
discussed here, under which the two poles in Eq. ([6]) are far apart from each other, and the eikonal approximation holds. 
We stress that to justify the eikonal approximation, fc+ needs not to be as small as kyr- Therefore, ordinary soft gluons 
exchanged between the active partons of, say, Hi and H2, can be still factorized using the eikonal approximation. In 
this case the parton, after emitting an ordinary soft gluon of the momentum Z^ ~ (A, A, A), carries the momentum 
fci — Z. We then have the hierarchy fc+Z _ 3> \\tix — 1t| 2 for fc^ 3> kiT,lr, which justifies the eikonal approximation 
l/(fci-Z) 2 «l/(-2fci-Z). 
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One may point out that the Wilson line for a TMD parton density should be rotated away from the light cone 
in order to regularize the light-cone singularity [25|, [2(| . Below we examine the effect of replacing the Wilson line 
direction = (0,1, Ot) by n M = {n + ,n~ ,0t) with n 2 < 0. The three denominators in Eq. ([5]), containing the 



terms 2{l+ - k+)l~ , 2(1+ + pf - kf)l' 
> 1+ > -(p+ -k+): 



and 2l + l , lead to three poles in different / half planes for the range 



r = fcr 



fei H 
2F 



•■IT I 



2(Z+ 



Ur-k 



IT | 



2(l++pt-kt) 



it, 



(10) 



(11) 



(12) 



Solving Eqs. (TTU1) . (TTTI) . and (|12l) together with the constraint n + Z + n l + = from the ^-function 5(n ■ Z), we obtain 



Z+ = I 
2 

I- = i 
2 

Z+ = - 



ki Pi k 1 



k++k- — -2 — 1 T 
n J n 



C1T 




Pi - k? - k l 



-2 — |l T -k 1T | 2 

n 



2n 



-h 



(13) 



(14) 



(15) 



respectively. It is then clear that the choice n 2 < prevents complex solutions. We have checked that no solution of l + 
exists in the range k^ > l + > 0, for which the three poles of l~ are also located in different half planes. Equations (TTU)) 
and (fl3|) imply the pole l~ ~ A 2 /fc^" + ie, if n + and n~ are of the same order of magnitude. Equations (|11[) and (fT4|) 
[ ([12"]) and (p"5j) ] imply the pole l~ ~ A 2 /E — it (l~ ~ A + ie). Hence, the contour of Z~ can be deformed, such that 
l~ remains of 0(A 2 /k^ ) in the contour integration. That is, the eikonal approximation for the spectator line of Hi 
is justified even under the variation of the Wilson line direction. 

We sum the residual soft corrections to all orders in the impact parameter space. The factors l/|kiT — Iy| 
(fe + fei) • (ki + k^)/(—2kfk^ — |kiT — k3T — 1t| 2 ) are absorbed into the LO TMD parton density fty and the 

LO parton-level differential cross section daf^^ k+l , respectively. The contribution from Figs, [lja) and [Ub) is then 
factorized into the convolution 



and 



T, 



(i) 



.9 



(2^ 



d 2 h 



rf^ ( + ) 7 _ >fe+; (kiT 



h 



(16) 



c 3T ,ki T - k 3 T)f^l(k 1T - lr,ki T ) 
is obvious. The argument kiT — It (kir) in / 7] labels the parton transverse 



where the soft divergence from It 

momentum in the hadron Hi before (after) the final-state cut. The two arguments in da 

active partons participate the hard scattering actually: one parton corresponds to the valence scalar of Hi, 
another to the Glauber gluon. Equation (|16p is rewritten, in the impact parameter space, as 



(0) 
i+j—tk+l 



indicate that two 
and 



T 



(i) 



d 2 b l d 2 b r d 2 b' l d 2 b' r [-iS(b l )]da 



(b I -b{ > b r -b^(b{ > b;) 



xe iki T -(b ! -b r ) e -'ik 3T -(b 1 -b r -bJ+b^ 



with the one-loop soft factor 



5(b) 



d 2 l 



T e -il T .h 



(17) 



(18) 



The factor —i has been made explicit, so that the soft factor 5 defined in Eq. (|T8|) is real. 

Here we give more explanation to the notations employed in Eq. (fT7| . The variables h\ and b' r denote the transverse 
coordinates of the partons coming out of Hi, namely, the lower ends of the virtual gluons, before and after the final- 
state cut, respectively. The definition of fan is given by the standard matrix element of the nonlocal operator 



fiMK) = I ^e- ix ^y'(Hi\4(y-,b' r )W-(y-,bUoo^W-(0, b{;oo)&(0,bD|iTi>, 



(19) 
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where the parton momentum fraction x\ has been omitted in Eq. (|17l) . The factor W- denotes the Wilson line 
operator 



W-(y , b;oo)=Pexp 



-19 



d\u- ■ A(y + Au_) 



(20) 



with the coordinate y = (0,y ,b). It should be understood that the two Wilson lines W-(y ,b' r y and W_(0,b|) 
are connected by a link at infinity [27| . The translational invariance of fai demands that it depends only on the 
difference b' r — b[ . The above notations correspond to the choice of the triple-scalar vertex in Hi as the origin of the 
transverse coordinates. Choices of other vertices as the origin are certainly allowed. Accordingly, b; and b r denote 
the transverse coordinates of the upper ends of the virtual gluons before and after the final-state cut, respectively. 
The factor S"(b/) in Eq. (|17p then describes the gluon propagation in the transverse plane from the triple-scalar vertex 
to the upper end of the virtual gluon before the final-state cut. 

If the soft gluons appear on the right-hand side of the final-state cut, we have 



T. 



(i) 



xe ik 1T .(b,-b r ) e 



-ik3 T -(bi-b r -h' t +h' r ) 



(21) 



At next-to-next-to-leading order, we have two radiative gluons either on each side or on the same side of the final-state 
cut. The former case leads to 



T, 



(2) 



LR 



d 2 6, d 2 M 2 MX H S (b z )] d<7, (0) 
Xe iki T .(b,-b r ) e -ik 3T .(b,-b r -b;+b;) _ _ , 



(22) 



For the latter with the two radiative gluons on the left-hand side of the final-state cut, we derive the factorization of 
the residual soft divergence 



T, 



(2) 



LL 



9 



(2tt)' 



'it 



'2T , (0) 



i2T 



£3Ti K 1T 



<-3T) 



(23) 



which can be rewritten as 



T, 



(2) 



LL 



d^Xd^dX^-iSib^dal^^ht - b;,b r - b'jftjKbiX) 



(0)/ 



ikiT-(bi-b r ) -ikaT-Cbi-b^-bJ+b;.) 



(24) 



To achieve at Eq. (f2U| from the summation of all the relevant two- loop diagrams, the following relation from the 
eikonal approximation for the spectator line of H\ must be employed 



1 



1 



1 



1 



1 



1 



+ it (l x + l 2 ) + if l 2 + if (li + l 2 ) + it li + ie 1% + «e 



(25) 



On the left-hand side of Eq. ([25|) . the second term is obtained by exchanging the gluon carrying the momentum l\ 
and the gluon carrying li in the first term. The above relation will not hold, if the corresponding transverse loop 
momenta are retained in the denominators. That is, the all-order factorization of the soft gluons from the process is 
not possible without the eikonal approximation. 

Viewing Eqs. ([23| and (|24|) . it is clear why the all-order summation of the residual soft divergences can be facilitated 
in the b convolution, instead of in the kx convolution. Applying the above procedure to higher loops, the diagrams 
with gluons being emitted by the spectator of the hadron Hi and attaching to the active partons of Hi and H4 are 
summed into 



dXdXifbidXe-^^dai+^k+tfa - b\,b r - b;)e jS ^/, ;/1 (b;,b;; 



x ikiT-(bi — b r ) e -ik3 r -(bi-b r 



- b ;+ b ;) . . . 



(26) 
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The collinear gluon exchanges of the type in Fig.QJc) can be factorized in the standard way, so the definitions of both 
d&i+j^-k+l an d fi/i have been extended to all orders. When the final-state hadron pair carries a net large transverse 
momentum, the k\x dependence in the fragmentation functions are negligible. Integrating over fax, the 5- function 
6 (hi — h r ) renders the soft factor vanish, and Eq. ((26)) reduces to a formula in the collinear factorization [28j . 
The operator definition of the soft factor is given by 

e -iS(b) = (o|T^_(o,b;-c») t W_(0,b;cx))W'+(0,0;c»)W+(0,0;-cx)) t |0) ) (27) 

where the factor W+ denotes another Wilson line operator 



W + (y + , b; oo) = P exp 



ig / d\u + ■ A(y + \u + ) 



(28) 



with the coordinate y = (y + ,0,h) and the dimensionless vector u+ = (1,0,0). The net effect of 
W-(0, b; — oo)^W_(0, b; oo) and W+(0, 0; oo)W+(0, 0; — ooy demands the vanishing of the components l + and l~ 
of the loop momentum, respectively. A Glauber gluon is then off-shell by l T as indicated in Eq. (fT8|) . It can be shown, 
by expanding the Wilson line operators order by order, that Eq. (|27p reproduces the Feynman rules obtained here for 
the soft factor. 

To derive the differential cross section, we integrate the momentum conservations 8(k^ — k% — k^), S(k^ — k% — k±), 
and 8 2 {k.iT + k 2 T — k 3 y — over , k^ and fax, respectively, and the on-shell conditions S(k^) and <5(fcf ) over 
/c 3 and k®, respectively. The longitudinal parton momenta k^ and k^ are then related to k3 and k4. Substituting 
kiT = k 3 T + k 4 T — k 2 T into the Fourier factor e lk iT-(b,-b,.) ^ an( j integrating over k 2 x, fj/2^2r) is transformed into 
the impact parameter space. At last, we arrive at the factorization formula modified by the soft factor associated 
with the hadron H\ 

E 3 E 4 „ d ° = V / f k , ^ d% l dXd 2 b' l d 2 b' r e- iS ^ ( B i+j ^ k+l (h l - bj,b r - b;)e'^ b ') 
cPp 3 cPp4 ^ J |k 3 | |k 4 | 

x/ t/1 (b;,b;)/ j/2 (b ; - b r )d fc/3 (k 3 )d i/ 4(k4)e ik3T - (b; - b; ' ) e 4k4T - (b! - b '- ) , (29) 

where the dependence on k 3 and k^ in d&i + j^k+i, fi/i and fj/2 is not shown explicitly. The residual soft divergences 
associated with the hadron H2 can be analyzed similarly, which are not discussed in this work. 

Note that the soft factor in Eq. (|27[) does not carry the flavor indices i, j, k, and I, and is independent of the 
species of hadrons involved in the collision. This universality makes possible experimental constraints on its behavior 
from some processes (e.g., HH — > tttt + X), and predictions from the kx factorization for other processes (e.g., 
HH — > KK + X). Equation (|29|) can also be employed to study the effect from the factorization breakdown: one can 
compare the numerical difference between the predictions from Eq. |29|) and from the formula without the soft factor. 
The soft factor derived above differs from that in the simple Drell-Yan process for which the kx factorization 
has been justified [l9|, |2(J. That is, the infrared divergences studied in [l(| arise from the ordinary (not Glauber) 
soft region, and can always be collected by means of the eikonalzation. The appearance of the soft factor in [l(| 
is attributed to the incomplete infrared cancellation between virtual corrections, where loop momenta do not flow 
through hard scattering, and real corrections, where loop momenta do. 

We next discuss the failure of the kx factorization for the transverse single-spin asymmetry (SSA) in hadron-hadron 
collision [29l - l33j with H\ being the transversely polarized hadron. For the SSA, the parton transverse momentum 
must be taken into account, and the imaginary part of the polarized TMD parton density contributes. Adopting a 
similar model field theory [3J], which contains additional fermion fields, the kx factorization for the SSA was also 
shown to fail [7J. The mechanism is identical to that in the unpolarized hadron hadroproduction. This is the reason 
why the definitions of the TMDparton densities in the SSA and unpolarized cases were modified by including the 
same additional Wilson links in [8j. The sum of these additional Wilson links leads to the (5-function in Eq. ©, which 
breaks the kx factorization at one loop. 

The loop integral associated with the sum of Figs, [lja) and [lib) f° r the SSA, where the solid lines represent 
fcrmions, is written as 



d 4 l (k 2 + fc 4 ) ■ (fci + k 3 )6(l+) 

(2tt) 4 (fci - lf( Pl -ki + /) 2 / 2 (-2fc+fc 3 - - |k 1T - k 3T - 1 T |2) ' 



K^Tr[(A+m ff ) 75 /(A- M+ /)l + (h- fa)]- (30) 

The spin vector s is chosen in the transverse direction, and tub denotes the mass of the hadron H\ . The trace in the 
above expression gives, in the small k\ region, 

^Tr[(h + m H ) l5 Mfa- fa+ t)l + (i>i- fa)\ ~ 2im H e a0 s a lPp+, (31) 
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where the e tensor obeys £12 = 1, and picks up the transverse components. Working out the integrations over / + 
and then over l~ in Eq. (|3"U)) . the denominator l 2 (ki — I) 2 becomes I t \\t — ^it\ 2 - Due to the existence of 1? in Eq. pip , 
a residual soft divergence arises from the region with lx — > k%x in this case. 

Following the same reasoning, we may be able to eikonalize the spectator quark line of Hi in the Glauber region 

1 1 1 1 - (32) 



(pi — k% + I) 2 + ie pi ■ I + it I + it 

where the numerator 2pf comes from Eq. pip . However, the rest of the procedure for factorizing the residual soft 
divergence in the unpolarized hadron hadroproduction does not apply: the propagator 1/|1t — kir| 2 can not be 
absorbed into the LO polarized TMD parton density, which vanishes with the fermion trace 

l -Tr[{ A + m H ) l5 Mfa- h)\ = 0. (33) 



It implies that the factorization of the residual soft divergence in Eq. (|30|) into the convolution of the one-loop soft 
factor, the LO polarized TMD parton density, and the LO hard kernel is not possible even at low px in this model 
field theory. In QCD, the corresponding TMD parton density, namely, the Sivers function, also vanishes at LO in 
perturbation theory. Therefore, the universality of the Sivers function can not be restored by factorizing out a soft 
factor. 

Our observation obtained in this work has wide applications. For example, the same residual soft divergences from 
the Glauber region have been identified in the color-suppressed tree amplitudes [35j j of two-body nonleptonic B meson 
decays in the perturbative QCD (PQCD) approach, which is based on the kx factorization [36l - l40l ] . A soft factor 
has been introduced into the PQCD factorization formulas, which enhances the color-suppressed tree amplitudes 
significantly, such that the known irn and nK puzzles were resolved [35j |. 

In a forthcoming paper we shall discuss the kx factorization of the residual soft divergences from low-px hadron 
hadroproduction in real QCD, where the spectator lines in Fig. Q] are replaced by infinitely many rung gluons. A 
difference arises from additional color degrees of freedom of quarks and gluons. Equation p3 is then a consequence 
of the summation over the attachments of the radiative gluon to the active parton lines in H2 and H4 and to the 
hard gluon line. To employ Eq. ([25]) in QCD, the diagrams with triple gluon vertices should be included too. We 
shall demonstrate, with the help of the eikonalization and the Ward identity, that the breakdown of the universality 
of a TMD parton density is suppressed in the large N c limit, N c being the number of colors, after summing over the 
attachments to all rung gluons. The kx factorization is also applicable to the W boson plus jet production and the 
direct photon production, for which the momentum ki (£3) is carried by a partonic gluon (gauge bosons). In these 
processes the (^-function leading to the residual soft divergence appears after summing over the attachments to the 
partonic gluon, the quark carrying the momentum fej, and the virtual quark. 
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